Genome-wide detection of selective signals for fecundity traits in goats (Capra hircus).

Fecundity in livestock is an economically important complex quantitative trait that is influenced by both genetics and the environment. However, the underlying genetic mechanism of reproductive performance in goats has not been well investigated. To investigate the genomic basis of fecundity in goats, genomic sequencing data of the Jining grey goat (a high prolificacy breed in China) were collected, as well as data for other commonly available goat breeds, and a mass of genomic variants were generated after variation calling. We screened the Jining grey goat (20 individuals) using a selective sweep with the Asian wild goat population (5 individuals), and potential candidate genes were proposed, such as STIM1, ESR1, LRRC14B and SLC9A3. Among, STIM1 is a most promising one associated with high reproductive capacity. When compared to Chinese domestic goats with low fecundity (17 individuals), the genes including MLLT10, SPIRE2, TCF25, ZNF276 and FANCA were screened, and the SPIRE2 gene was thought to be associated with fecundity traits. Meanwhile, the functional enrichment of these candidate genes revealed that they were involved in biological processes of mammary gland morphogenesis, uterus development, gastrulation, mesoderm morphogenesis and formation, and blood vessel development, which might undergo natural or artificial selection during reproductive trait formation in goats. Thus, our findings could enrich the genetic basis of reproductive trait selection during goat domestication, which may serve to improve goat breeding practices.

Genome-wide association study identifies quantitative trait loci affecting cattle temperament.

Cattle temperament is an interesting trait due to its correlation with production efficiency, labor safety, and animal welfare. To date, however, its genetic basis is not clearly understood. Here, we performed a genome-wide association study for a series of temperament traits in cattle, assessed with via open field and novel object tests, using autosomal single nucleotide polymorphisms (SNPs) derived from the whole-genome sequence. We identified 37 and 29 genome-wide significant loci in the open field and novel object tests, respectively. Gene set analysis revealed the most significant pathway was the neuroactive ligand-receptor interaction pathway, which may be essential for emotional control in cattle. Analysis of the expression levels of 18 tissue-specific genes based on transcriptomic data showed enrichment in the brain, with some candidate genes involved in psychiatric and neurodegenerative diseases in humans. Based on principal component analysis, the first principal component explained the largest variance in the open field and novel object test data, and the most significant loci were assigned to SORCS3 and SESTD1, respectively. Our findings should help facilitate cattle breeding for sound temperament by pyramiding favorable alleles to further improve cattle production.

Mitochondrial genomes from modern and ancient Turano-Mongolian cattle reveal an ancient diversity of taurine maternal lineages in East Asia.

Turano-Mongolian cattle are a group of taurine cattle from Northern and Eastern Asia with distinct morphological traits, which are known for their ability to tolerate harsh environments, such as the Asian steppe and the Tibetan plateau. Through the analysis of 170 mitogenomes from ten modern breeds, two sub-lineages within T3 (T3119 and T3055) were identified as specific of Turano-Mongolian cattle. These two T3 sub-lineages, together with the previously identified T4, were also present in six Neolithic samples, dated to ~3900 years BP, which might represent the earliest domestic taurine stocks from Southwest Asia. The rare haplogroup Q, found in three Tibetan cattle, testifies for the legacy of ancient migrations from Southwest Asia and suggests that the isolated Tibetan Plateau preserved unique prehistoric genetic resources. These findings confirm the geographic substructure of Turano-Mongolian cattle breeds, which have been shaped by ancient migrations and geographic barriers.

Single-cell transcriptome landscape of ovarian cells during primordial follicle assembly in mice.

Primordial follicle assembly in the mouse occurs during perinatal ages and largely determines the ovarian reserve that will be available to support the reproductive life span. The development of primordial follicles is controlled by a complex network of interactions between oocytes and ovarian somatic cells that remain poorly understood. In the present research, using single-cell RNA sequencing performed over a time series on murine ovaries, coupled with several bioinformatics analyses, the complete dynamic genetic programs of germ and granulosa cells from E16.5 to postnatal day (PD) 3 were reported. Along with confirming the previously reported expression of genes by germ cells and granulosa cells, our analyses identified 5 distinct cell clusters associated with germ cells and 6 with granulosa cells. Consequently, several new genes expressed at significant levels at each investigated stage were assigned. By building single-cell pseudotemporal trajectories, 3 states and 1 branch point of fate transition for the germ cells were revealed, as well as for the granulosa cells. Moreover, Gene Ontology (GO) term enrichment enabled identification of the biological process most represented in germ cells and granulosa cells or common to both cell types at each specific stage, and the interactions of germ cells and granulosa cells basing on known and novel pathway were presented. Finally, by using single-cell regulatory network inference and clustering (SCENIC) algorithm, we were able to establish a network of regulons that can be postulated as likely candidates for sustaining germ cell-specific transcription programs throughout the period of investigation. Above all, this study provides the whole transcriptome landscape of ovarian cells and unearths new insights during primordial follicle assembly in mice.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences.

BACKGROUND: Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal. RESULTS: In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits. CONCLUSIONS: The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

Genomic Signatures of Selection Associated With Litter Size Trait in Jining Gray Goat.

Litter size (LS), an important economic trait in livestock, is so complicate that involves many aspects of reproduction, the underlying mechanism of which particularly in goat has always been scanty. To uncover the genetic basis of LS, the genomic sequence of Jining Gray goat groups (one famous breed for high prolificacy in China) with LS 1, 2, and 3 for firstborn was analyzed, obtaining 563.67 Gb sequence data and a total of 31,864,651 high-quality single nucleotide polymorphisms loci were identified. Particularly, the increased heterozygosity in higher LS groups, and large continuous homozygous segments associated with lower LS group had been uncovered. Through an integrated analysis of three popular methods for detecting selective sweeps (Fst, nucleotide diversity, and Tajima's D statistic), 111 selected regions and 42 genes associated with LS were scanned genome wide. The candidate genes with highest selective signatures included KIT, KCNH7, and KMT2E in LS2 and PAK1, PRKAA1, and SMAD9 in LS3 group, respectively. Meanwhile, functional terms of programmed cell death involved in cell development and regulation of insulin receptor signaling pathway were mostly enriched with 42 candidate genes, which also included reproduction related terms of steroid metabolic process and cellular response to hormone stimulus. In conclusion, our study identified novel candidate genes involving in regulation of LS in goat, which expand our understanding of genetic fundament of reproductive ability, and the novel insights regarding to LS would be potentially applied to improve reproductive performance.

Determination of genetic effects of SERPINA3 on important growth traits in beef cattle.

Serine protease inhibitor protein 3 (serpin peptidase inhibitor, clade A, member 3, SERPINA3) is a member of the serpin superfamily, probably related to the yield and quality of muscle. This study focuses on the relationship between SERPINA3 gene polymorphism and growth traits in beef cattle. The study first uses sequencing pooled DNA samples (Pool-Seq), PCR-RFLP and Tetra-primer ARMS-PCR techniques to determine the genetic polymorphisms of SERPINA3 in 765 beef cattle. Then, the polymorphic loci were correlated with the growth characters of cattle. Five SNPs (SNP1:A-648G, SNP2:T6496A, SNP3:G2495A, SNP4:T2595A, SNP5:A2615G) were found, located in the promoter, introns 5 and SNP 3, 4, 5 were in exons 2, respectively. The observed He was from 0.44 to 0.5, Ne were approaching 2 (1.78 to 2.00). The maximum and minimum PIC (polymorphism information content) values were 0.37 and 0.34, respectively. The association analysis results showed that the SNPs had a significant height in the chest girth and body length. (p < 0.05 or p < 0.01). This will provide important information for the rapid breeding of Chinese yellow cattle and the establishment of a molecular genetic marker database.

Detection of InDel variations within seven candidate genes and their associations with phenotypic traits in three cattle breeds.

The Xinjiang brown cattle, Red steppe cattle, and Yunling cattle are indigenous cultivated cattle breeds in Chinese frontier provinces, and they produce high-grade beef and milk products, however, their genetic diversity in many important genes related to excellent meat and milk production is still unknown. Our previous studies have found that several candidate genes (e.g., SREBP1c and PAX7) were associated with bovine economically important phenotypic traits, but none has been reported in the above-mentioned three cattle breeds. Since the InDel (insertion/deletion) marker becomes a useful tool applied in the animal molecular breeding, herein, we firstly found that the InDel variations of seven candidate genes in these cattle. Results showed that the genotypic and allelic distributions of these seven genes were remarkably different among these three cattle (p < 0.05 or p < 0.01). Furthermore, the InDel variations of SREBP1c and PAX7 genes were significantly associated with eight phenotypic traits in Xinjiang brown cattle (p < 0.05 or p < 0.01), respectively, suggesting that they can become the useful DNA markers.

Integrative analysis of APOL3 gene CNV for adult cattle stature.

Copy number variations (CNVs) have been identified as another important structural variation of genome. In recent years, a large amount of CNVRs have been identified in humans and animals. However, association and dosage effects studies of CNVs are very limited. Apolipoprotein L3 (APOL3) gene plays a central role in modulating gene transcription and is located within a CNVR that encompasses quantitative trait locis (QTLs) for economic traits like meat quality. Herein, we analyzed the CNV polymorphism of APOL3 in 421 individuals from five distinct cattle breeds, and then correlated their genotypes with growth traits. Association analysis revealed that the APOL3 CNV was significantly associated with hip height and cannon circumference of Xianan (XN) cattle (P < .01), and visibly associated with body slanting length and hucklebone width of Pinan (PN) cattle (P < .05). Overall, the data provide evidence for the functional role of APOL3 CNV and a basis for future applications in cattle breeding.

Copy number variation of MYLK4 gene and its growth traits of Capra hircus (goat).

Copy number variation (CNV) is a form of genetic variation caused by genome rearrangement, with abnormal fragments ranging from 50 bp to Mb. And, CNV is closely related to disease, growth and reproductive shape of livestock. As a member of myosin light chain kinase (MYLK) family with serine/threonine specificity, MYLK4 belongs to an enzyme encoded by MYLK4 gene. Although MYLK4 is a recognized kinase, its function has yet to be revealed in subsequent studies. This study aims to analyze CNV and genetic effects of MYLK4 gene in goats. We used qPCR to detect CNV of MYLK4 gene in African Nubian goat (n = 32), Guizhou black goat (n = 196) and Guizhou white goat (n = 95), respectively, and correlated CNV data of MYLK4 gene with goat growth traits in Chinese goats. The results showed that the effect of MYLK4 gene CNV on body weight, body length and body height of goats had significantly different (p < 0.05, Q < 0.05), in which CNV showed better growth traits in type of deletion. Therefore, CNV of MYLK4 gene can be used as a molecular marker for assisted selection of goat growth traits, which provides a theoretical basis for the genetic improvement of goat breeds in China.

Exploring genotype-phenotype relationships of the CRABP2 gene on growth traits in beef cattle.

Cellular retinoic acid binding protein 2 (CRABP2) is essential to myoblast differentiation. However, there was no report about the function of CRABP2 gene in cattle. This study explored the association of CRABP2 gene polymorphisms with growth traits in cattle breeds by several methods, such as DNA sequencing, PCR, PCR-RFLP and forced PCR-RFLP. Two sequence variants were determined. There were 621 individuals in six cattle breeds from China for the experiment, and three breeds were used to test validation of polymorphisms and extent of linkage disequilibrium (LD). The results showed that both SNPs (SNP1, g.2458 G > T, SNP2, g.3878 G > A) were in intron1. Two SNPs were in low linkage disequilibrium. Association analysis suggested that SNP1 had the significant difference on growth traits with body height, height at hip cross and body slanting length (P < .05), while SNP2 showed a significant difference in growth traits with body height, height at hip cross and body slanting length(P < .05). The results of this investigation displayed that the CRABP2 gene is an available candidate gene and may be used for breed selection and conservation.

Copy number variation (CNV) in the IGF1R gene across four cattle breeds and its association with economic traits.

The insulin-like growth factor 1 receptor (IGF1R) plays a vital role in immunomodulation and muscle and bone growth. The copy number variation (CNV) is believed to the reason for many complex phenotypic variations. In this paper, we statistically analyzed the copy number and the expression profiling in different tissue types of the IGF1R gene using the 422 samples from four Chinese beef cattle breeds, and the mRNA of IGF1R was widely expressed in nine tissue types of adult cattle (heart, liver, kidney, muscle, fat, stomach, spleen, lung and testis). Results of CNV and growth traits indicated that the IGF1R CNV was significantly associated with body weight and body height of Jinnan (JN) cattle and was significantly associated with body height and hucklebone width of Qinchuan (QC) cattle, making IGF1R CNV a promising molecular marker to improve meat production in beef cattle breeding. Bioinformatics predictions show that the CNV region is highly similar to the human genome, and there are a large number of transcription factors, DNase I hypersensitive sites, and high levels of histone acetylation, suggesting that this region may play a role in transcriptional regulation, providing directions for further study of the role of bovine CNV and economic traits.

Copy Number Variation of the SHE Gene in Sheep and Its Association with Economic Traits.

Copy number variation (CNV) caused by gene rearrangement is an important part of genomic structural variation. We found that the copy number variation region of the Src Homology 2 Domain Containing E (SHE) gene correlates with a quantitative trait locus of sheep related to milk fat percentage and bone density. The aim of our study was to detect the copy number variation of the SHE gene in four sheep breeds and to conduct a correlation analysis with economic traits, hoping to provide some reference for sheep breeding. In this study, we examined 750 sheep from four Chinese breeds: Chaka sheep (CKS), Hu sheep (HS), Large Tail Han sheep (LTHS) and Small Tail Han sheep (STHS). We used qPCR to evaluate the copy number of the SHE gene, and then used general linear models to analyze the associations between CNV and economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS individuals. Association analyses showed that gain and normal copy number types were correlated to body length, circumference of cannon bone, heart girth, chest width and high at the cross in CKS, HS and STHS (p < 0.05), but this association was not observed for LTHS. Chi-square values (χ2) found prominent differences in CNV distribution among the studied breeds. Overall, the CNV of the SHE gene may be an important consideration for sheep molecular breeding.

Genetic Variants in ADD1 Gene and their Associations with Growth Traits in Cattle.

The α-adducin (ADD1) is a subunit of adducin which is a cytoskeleton heterodimeric protein. Adducin participates in oocytes chromosome meiosis of mice, prompting adducin has an effect on embryonic development. Adducin gene mutation has significantly functional change. So the present study was to identify and characterize polymorphisms within the coding region of the bovine ADD1 gene among different cattle breeds. Here, 11 novel single nucleotide polymorphisms (SNPs 1-11) were identified by DNA sequencing and polymerase chain reaction-single stranded conformational polymorphism, there were one synonymous mutation in exon 1 (SNP1); four missense mutations in exons 4, 7, and 8 (SNPs 3-6); and six mutations in introns 4, 12, 13, and 14 (SNPs 2, 7-10). The statistical analyses indicated that the some SNPs are associated with the growth traits (body length, body height, chest circumference, and hucklebone width) in Chinese Jiaxian cattle population. Our results provide evidence that polymorphisms in the ADD1 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.

Novel copy number variation of the KLF3 gene is associated with growth traits in beef cattle.

Copy number variation (CNV) related to complex traits, such as disease and quantitative phenotype, is considered an important and wealthy source of genetic and phenotypic diversity. It suggests that the copy number variation of function gene maybe leads to the phenotypic changes. Kupple like factor 3 (KLF3) gene is a vital transcription factor associated with the growth and development of muscle and adipose tissue. It has been mapped in a CNV region by animal genome re-sequencing. In this study, we detected the distribution diversity of KLF3 gene copy numbers in six Chinese cattle breeds (QC, NY, XN, PN, QDM and JX) and associated the phenotypic traits with it. Then, we analyzed the KLF3 gene transcription expression level in different tissues of Jiaxian (JX) cattle. Furthermore, we detected mRNA expression level of muscle and fat tissues of Jiaxian cattle (JX), Angus × Jiaxian (AJ). The results showed that the copy number in CNV loss was more frequent in QC than others. And we revealed a positive effect of KLF3 CNV on growth traits, such as body mass and heart girth (P < 0.05). In a word, we ascertained the significance between CNVs of KLF3 gene and growth traits in different cattle breeds, and our data indicates that the CNVs of KLF3 gene may as a marker for the future molecular breeding of Chinese beef cattle.

Genetic Variant of MYLK4 Gene and its Association with Growth Traits in Chinese Cattle.

As a member of MYLK family, MYLK4 gene may play a vital role in muscle development. In this study, one novel single-nucleotide polymorphism (SNP) was identified the bovine MYLK4 by sequencing pooled DNA samples (pool-Seq) and forced polymerase chain reaction-restriction fragment length polymorphism (forced PCR-RFLP) methods. Overall, we reported one mutation (SNP1) in the intron 10 region within the bovine MYLK4 gene in 559 individuals representing five main cattle breeds from China (Nanyang, NY; Qinchuan; Jiaxian, JX; Pinan cattle; and Caidamu cattle, CDM). Genotype AA and allele A were predominant in the QC, PN, and XN populations. Association analysis with growth traits in the QC breed showed that the animals with genotype GG had significantly greater chest breadth and hip width (P < 0.05). Meanwhile, the genotype GG was strongly associated with withers height and body length than those with genotype AA (P < 0.01 or P < 0.05) at 12 months in the NY breed. These statistical results exhibited that the MYLK4 gene might be a potential candidate gene to improve cattle's growth traits, and the SNP could be used as molecular markers in early marker-assisted selection (MAS) in beef cattle breeding program.

Distribution and association study in copy number variation of KCNJ12 gene across four Chinese cattle populations.

Copy number variation is a large genome variation which usually happens in the noncoding-region, and it may occur at the locus associated with the functional gene to further influence the phenotype. Potassium inwardly-rectifying channel, subfamily J 12 (KCNJ12) gene expressed widely in cardiomyocytes and neurons, plays an important role in tumor therapy and muscle movement regulation. In this study, we detected the distribution of CNVs for KCNJ12 gene in 404 individuals belonging to four Chinese cattle breeds (NY, JX, JA and GF). We also investigated the KCNJ12 gene expression in different tissues of JX cattle. Additionally, we examined the association of two CNV regions (CNV1: 1,600 bp, intron 1; CNV2: 4,800 bp, intergenic) with growth traits. The statistical analyses indicated that the CNV1 is associated with the body length, rump length and weight in JX cattle population (P < 0.05); and there has a significant association with the body length, chest circumference, and body weight in GF cattle (P < 0.05).The CNV2 had a significant effect on the body length and body weight in JX cattle (P < 0.05); the body length, chest circumference, rump length and body weight in GF cattle (P < 0.01 or P < 0.05). The copy numbers of KCNJ12 gene presented the negative correlations with the transcript level of gene in skeletal muscles (P < 0.05). Our results provide evidence that CNV1 and CNV 2 in KCNJ12 are associated with growth traits in two cattle populations and may be used as candidates for marker-assisted selection and breeding management in cattle.

A novel SNP of PLAG1 gene and its association with growth traits in Chinese cattle.

Genetic polymorphism has great influences on the improvement of cattle traits. The polymorphism of the same gene family will greatly change the growth traits of cattle, such as the pleomorphic adenoma gene (PLAG) family. Many studies have shown that the PLAG family proteins are the transcription regulators of nuclear protein, which mainly regulates the expression of many important genes in the body. In cattle, single nucleotide polymorphisms (SNPs) within or near the PLAG1 gene is associated with economic traits of height in cattle. Here we investigate a SNPs of bovine PLAG1 in 646 cattle from five breeds. We found three different genotypes by using cleaved amplification polymorphism sequence-tagged sites (CAPs), includes some significant differences in body height, chest circumference and other shapes (P < 0.05), also we found that the TT genotype had no advantage in body shape. These results indicate that the selection of PLAG1 gene could be used to ensure the breeding direction for growth traits of the beef cattle.

DNA methylation status of CRABP2 promoter down-regulates its expression.

As an important epigenetic modification DNA methylation is catalyzed by DNA methylation transferases (DNMTs) and occurs mainly in CpG islands. DNA methylation plays an important role in regulates gene expression, cell differentiation, genetic imprinting and tumor therapy. Retinoic acid-binding proteins (RAC) is vital for the absorption, transport, metabolism and maintenance of homeostasis of retinoic acid, which in turn regulates the differentiation and proliferation of cells by regulating the transcription of many target genes, therefore, these proteins influence differentiation and proliferation of adipocytes and muscle fibroblasts. Thus, cellular retinoic acid binding protein 2 (CRABP2) may be a candidate gene which affects beef quality, yield and fat deposition. The aim of this study was to evaluate the expression and the methylation pattern on the differentially methylated region (DMR) of the promoter of CRABP2. The DNA methylation pattern was tested by bisulfite sequencing polymerase chain reaction (BSP), the quantitative real-time PCR (qPCR) was used to analysis the expression of CRABP2 gene. The results showed that the DNA methylation level was higher in purebred cattle breed than that in hybrid cattle breeds which was negative correlation with the expression of the gen. These results indicate that the methylation status of the CRABP2 DMR can regulate mRNA expression. What's more, there are different methylation and expression patterns in different breeds and tissues which may influence the phenotype, and the results may be a useful parameter to investigate the function of CRABP2 in muscle and fat developmental in Chinese cattle.